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In this paper we present a new algorithm for a layout optimization problem: this concerns 
the placement of weighted polygons inside a circular container, the two objectives being to 
minimize imbalance of mass and to minimize the radius of the container. This problem carries 
real practical significance in industrial applications (such as the design of satellites), as well as 
being of significant theoretical interest. Previous work has dealt with circular or rectangular 
objects, but here we deal with the more realistic case where objects may be represented as 
polygons and the polygons are allowed to rotate. We present a solution based on simulated 
annealing and first test it on instances with known optima. Our results show that the algorithm 
obtains container radii that are close to optimal. We also compare our method with existing 
algorithms for the (special) rectangular case. Experimental results show that our approach 
out-performs these methods in terms of solution quality. 



1. Introduction 

The Layout Optimization Problem (LOP) concerns the physical placement of in- 
struments or pieces of equipment in a spacecraft or satellite. Because these objects 
have mass, the system is subject to additional constraints (beyond simple Cartesian 
packing) that affect our solution. The two main constraints that we handle in this 
paper are (1) the space occupied by a given collection of objects (envelopment), 
and (2) the non-equilibrium (i.e. imbalance) of the system. The rest of the paper is 
organized as follows: In Section 2 we first present a detailed description of the prob- 
lem, and describe previous related research. In Section 3 we describe our algorithm, 
and in Section 4 we give the results of numerical experiments. 

2. The Layout Optimization Problem in Satellites 

The Layout Optimization Problem (LOP) was proposed by Feng et al. (@) in 1999, 
and has significant implications for the cost and performance of devices such as 
satellites and spacecraft. It concerns the two dimensional physical placement of 
a collection of objects (instruments or other pieces of equipment) within a space- 
craft/satellite "cabinet", or container. The LOP is demonstrably NP-hard (@). Early 
work on this problem (fill . [lil. [Til. [l7l) almost always modeled objects as circles in or- 
der to simplify the packing process. However, in real-world applications, objects are 
generally rectangular or polygonal in shape, and modeling them as circles leads to 



t Email: M.Amos@mmu.ac.uk 

ISSN: 0305-215X print/ISSN 1029-0273 online 

© 200x Taylor & Francis 

DOI: 10.1080/0305215YYxxxxxxx 

http://www.informaworld.com 



September 29, 2008 16:22 Engineering Optimization submission 



expensive wastage of space. We have recently reported work on solving the rectan- 



gular case (114 ) . and here we report a new algorithm (based on a different approach) 
to solve the polygonal case. 

We now briefly introduce related work on the packing of irregular items. Dowsland 
et al. use a so-called "Bottom-Left" strategy to place poly gona l items in a bin 
(0), with items having fixed orientations. Poshyanonda et al. (fl3) combine genetic 
algorithms with artificial neural networks to obtain reasonable packing densities. 
In other related work, Bergin et al. study the packing of identical rectangles in 
an irregular container |l|, 0). Burke and Kendall have applied simulated annealing 
to translational polygon packing (i.e., without rotations) (0). Other authors have 
applied simulated annealing to solve the problem of rotational polygon packing on 
a continuous space (0, H) . 

Given the additional constraint that imbalance of mass must be minimised, it is 
difficult to see how these existing methods may be directly applied to the current 
problem. In what follows we describe a new nonlinear optimization model for the 
LOP, and then show how it may be solved using simulated annealing. 



2.1. Notation and Definitions 

Here we describe the formal optimization model, by first explaining our notation 
for the representation of polygons. We then show how to quantify relations between 
polygons (such as distance and degree of overlap), which are central to the problem 
of assessing the overall quality of a layout. 

2.1.1. Structure of a polygon 

Suppose there are k polygons (1,2,..., k) to be packed. The structure of a polygon 
includes both its shape and its mass. We use str(i) to denote the initial structure 
of a polygon, i: 



str(i) = (ni,m.i, (n,r 2 , ...,r n J, (0 1 ,0 2 , -,0m)) (1) 

where rrii is the mass of polygon i, and n, is the number of vertices in the graph 
representation of polygon i. The positions of the n« vertices are defined by two lists 
of polar coordinates. List (r\,T2, ■■■, r m) defines the Euclidean distance from each of 
the rtj vertices to the polygon's centre of mass, and list (0i,02, ■-,0n i ) defines the 
orientation of each of the n, vertices relative to the centre of mass. Figure Q] shows 
how to define the initial structure of a square with edge length 1; in Figure Q] (a), 
the shape's centre of mass is located at the shape centre, whereas in Figure [Hb), 
the centre of mass is located at one vertex. We define the point of reference of each 
polygon as its center of mass, in order to simplify the notation. 

2.1.2. Radius of a polygon 
The radius of polygon i is defined as the maximum of (rj , T'i . . . . , r Ui ) : 



r(i) = max{ri,r 2 ,...,r rii } (2) 



With the polygon's centre of mass at its own centre, the circle with radius r(i) 
defines the minimum-sized circle that can completely cover the polygon. 
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2.1.3. State of a polygon 

We use Cartesian coordinates to record the positions of the polygons, and set the 
center of the container (that is, the circle) as the original point. We use sta(i) to 
denote the state of a polygon i: 



sta(i) = (xi,yi,a>i) (3) 

where Xi,yi is the position of the centre of mass, and a defines a rotation angle. 
Then with str(i) and sta(i), we can draw a polygon, i, as in Figure [21 

2.1.4. Distance between two polygons 

The distance between two polygons is defined as the Euclidean distance between 
their centres of mass: 



dis(i,j) = yj (xi - xj) 2 + (yi - yj) 2 ) (4) 

2.1.5. Overlap between two polygons 

If two polygons do not overlap, this measure is zero. If two polygons i and j 
overlap at all, we measure this as: 
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Figure 3. The overlap function is not continuous 



ove(i,j) = max{0, r(i) +r(j) — dis(i,j)} (5) 
This measurement of overlap has certain characteristics: 

• In Equation j5j), r(i) and r(j) are two constants, therefore dis(i, j) can be directly 
obtained by the positions of the two polygons. 

• It is clear that ove(i,j) > 0. To satisfy the non- overlapping constraint, we should 
minimize ove(i,j) to zero. 

• ove(i,j) is not a continuous function of the positions. As shown in Figure [H 
when two squares with edge length 2 are adjacent on one side, their overlap is 
zero, but when the left square is moved a little to the right, the overlap "jumps" 
to 2^/2 -2. 

Ascertaining overlap between two polygons is not a difficult problem in compu- 
tational geometry or computer graphics. In this paper, we look at each edge of one 
polygon in turn; if it is intersected by any edge of another polygon, then an overlap 
exists; otherwise, if one polygon is contained within another, then clearly an overlap 
exists. So the ascertaining of overlap has complexity 0(mn), where m,n are the 
number of edges of the two polygons. 

2.1.6. State of a layout 
A layout X is defined as the combination of the states of k polygons: 



x = {xi,yi,Oii,X2,y 2 ,<X2, ...,x k ,yk,ak) (6) 

2.1.7. Radius of a layout 
If (X x ,X y ) denotes the position of the centre of mass of a layout X, then 



X x 



vn%x% 



rru 



,Xy 



2^i=i 



rrnyi 



(7) 



We define the radius r(X) of a layout X as the longest Euclidean distance from 
its centre of mass to any of the vertices of the polygons. Because of the imbalance 
constraint, we place the centre of mass of the layout at the container center, so 
r(X) defines the minimum-sized container. 

2.1.8. Overlap of a layout 
The overlap of a layout is the sum of all overlaps between its polygons: 
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(8) 



i=l j=\,j^i 



2.1.9. Problem definition 

From the definitions above, we obtain an unconstrained optimization problem: 



where Ai, A2 are two constants. Because the overlap function ove(X) is not con- 
tinuous, f(X) is not continuous. In general, the overlaps are extremely deleterious, 
so Ai should be set large enough to prevent their introduction. However, we note 
that the computation does not introduce overlaps when attempting to decrease the 
radius of the layout at the final stage of optimization, because of the discontinuous 
of the ove function. 



3. Simulated annealing algorithm 

Simulated annealing (SA) is a probabilistic meta-heuristic that is well-suited to 
global optimization problems (|10l ). Annealing refers to the process of heating then 
slowly cooling a material until it reaches a stable state. The heating enables the 
material to achieve higher internal energy states, while the slow cooling allows the 
material more opportunity to find an internal energy state lower than the initial 
state. SA models this process for the purposes of optimization. A point in the 
search space is regarded as a system state, and the objective function is regarded 
as the internal energy. Starting from an initial state, the system is perturbed at 
random, moving to a new state in the neighbourhood, and a change of energy A.E 
takes place. If AE <0, the new state is accepted (a downhill move), otherwise the 
new state is accepted with a probability exp( ^j? ) (an uphill move), where T is 
the temperature at that time and is Boltzmann's constant. When the system 
reaches equilibrium, T is decreased. When the temperature approaches zero, the 
probability of an "uphill" move becomes very small, and SA terminates. 

Let to denote the initial temperature, imax denote the maximum number of itera- 
tions, E(x) denote the energy function, and emax denote the "stopping" energy.The 
general pseudo-code for simulated annealing may be written as follows: 



Algorithm 1: Standard SA 

set the initial state x and initial temperature t = to, let i = 
while i < imax and E(x) > emax 

perturb x in its neighbourhood and get x' 

if E(x') < E(x) then x = x' 

else x = x' with probability exp( E ^ t E( ~ x ^ ) 

decrease the temperature 

i = i+l 
return x 



minimize f(X) = \\ove{X) + \2r{X) 



(9) 



The performance of SA may be affected by several parameters: the initial temper- 
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ature t , the maximum number of iterations, imax, the "stopping" energy, emax, 
the structure of the neighbourhood, and the schedule of cooling. For a given prob- 
lem, the values of these parameters should be carefully selected. 



3.1. SA for polygon packing 

3.1.1. Neighbourhood structure 

In each iteration of SA, we perturb one polygon, thus obtaining a new layout, and 
then decide whether to accept or reject the new layout by means of an evaluation. 
In the ith iteration, the (i mod k)th polygon will be perturbed. Given an initial 
radius i?o, which is large enough to contain the polygons, the neighbourhood for 
polygon j is defined as: 



Xj,yj G ( h 1.05) x R x randomi— 1, 1) (10) 

i — 2x imax 



a, e ( 1- 1.05) x 7r x randomi- -1, 1) (11) 

J i — 2 X imax 

Equations (fTOl) and (fTTj) show that, at the beginning of the algorithm's execution, 
the position of a polygon may vary by (— 0.55-Ro> 0.55-Ro), and its orientation may be 
perturbed by (— 0.557T, 0.557r). This neighbourhood is large, and and the polygon can 
thus "explore" more space. As the algorithm proceeds, the neighbourhood becomes 
increasingly smaller. At the end of the algorithm's execution, the neighbourhood 
shrinks to 0.05 times its original size, then SA chooses the best solution in the 
neighbourhood. 

3.1.2. Temperature decreasing 

We use a simple rule to decrease the temperature: every cmax iterations, we let 
t = d x t, where d < 1. 

3.1.3. Description of the algorithm 

The detailed SA algorithm for polygon packing is therefore described as follows: 

Algorithm 2: SA for the packing problem 

randomly generate an initial layout X. Let t = to, i = 
while i < imax and f(X) > emax 
let j = (i mod k) 

randomly select Xj , i/j , ctj by ([TO]) and (fTTl) and get new X' 
if f(X') < f(X) then X = X' 

else X = X' with probability exp ^^"^^ ^ ) 
if i mod cmax = cmax — 1 then t = d x t 
i = i + l 
return x 
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4. Numerical Results 

We are not aware of any standard library of benchmark instances for this partic- 
ular problem, although such libraries do exist for other related problems (|4j). We 
therefore take a two-stage approach to testing our algorithm; we first design six in- 
stances with known optima, against which we may initially validate our method. We 
note that these instances include both convex and nonconvex polygons. After es- 
tablishing the effectiveness of our SA algorithm, we then test our algorithm against 
other recently-described methods for rectangle packing (rectangles, of course, being 
members of the polygon class), using both existing instances from the literature 
and new, larger instances. 



4.1. Known Optima 

The instances with known optima are described in Table [U, with graphical repre- 
sentations given in Figured! 

For each instance, we use SA to try to find the optimal layout. The value of imax 
is set to 20000 x k, the value of cmax is set to 100 x k, and the initial temperature 
is set to 100. The constants Ai and A2 are set to 100 (to induce a large f(x) and 
adjust the probability of uphill movement). In each case, the algorithm is executed 
40 times. 

Values for the best radius found, rbest, mean radius, f, and variance v are pre- 
sented in Table (2) Representations of the best results obtained are given in Figure 
From Table [2] and Figure 02 we observe that the SA algorithm can find layouts 
that are very close to the the optimal configuration for instances 1, 2, and 3, where 
the number of polygons is relatively small and the overall structures are simple. 
The optimal radius for instance 1 was originally calculated as f y^ = 2.121, but 
the our results yielded a smaller value. This prompted a re-estimation, giving a 

new optimum of y 4^). In the first three instances, the errors to the optimum are 
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(1) (2) (3) 




(4) (5) (6) 

Figure 4. Instances with known optima. 



about rbas * = 2%. The algorithm performs less well on instance 4, with 12 

polygons in a relatively complex configuration. In this case, our algorithm cannot 
find the best configuration, and the error to the optimum is 15%. Instances 5 and 6 
feature noncovex polygons. Because of the complexity of the shapes, the algorithm 
is unable to solve these instances to optimality. The errors to the optimal radius 
are 11% and 4% for instances 5 and 6 respectively. 



Table 2. Numerical results for SA on four instances with known optima 
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4.2. Rectangular Instances 

We now test our algorithm on four instances of the LOP containing only rectangu- 
lar shapes. The first three instances were first described in (flil ). and the fourth in 
(0), 

where all four instances were used to benchmark three different approaches: 
a genetic algorithm (GA), particle swarm optimisation (PSO) and a hybrid com- 
paction algorithm followed by particle swarm local search (CA-PSLS). Depictions 
of each instance are depicted in Figure El and full descriptions are given in Table 
[3) Since both particle-based algorithms out-performed the GA, we do not consider 
this last method here. 

We run each algorithm 50 times on each instance, recording the best radius 
found, rbest, average radius, r, standard deviation of radii, r a and average run time 
in seconds, t. Each algorithm is coded in C, compiled with g++ 4.1.0 , and run 
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■1 



Figure 6. Four instances of the LOP using rectangles 

under SUSE Linux 10.1 (kernel version 2.6.16.54-0.2.5-smp) on a computer with 
dual Intel Harpertown E5462 2.80GHz processors, 4GB of RAM and an 80GB hard 
drive. 

The three algorithms (SA, CA-PSLS and PSO) each run over a number of it- 
erations, which is dictated by the value of the constant CYCLE. In this set of 
experiments, we set CYCLE=3000 for each algorithm. The SA parameter values 
for cmax, initial temperature, Ai, and A2 are set as before, and the imax values set 
as to 100000, 120000, 108000 and 100000 for instances 1-4 respectively. The results 
obtained are depicted in Table HI 

On the first three (small) instances, both particle-based algorithms slightly out- 
perform the SA method in terms of solution quality; on average, by 10%. However, 
this comes at a significant cost disadvantage in terms of run time; over the first 
three instances, the particle-based methods require four times the execution time 
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Table 3. Four instances from the literature 

Instance k J?q Structure 

1 5 20 str(l) = (4, 12, (v"25, V25, \/25, V25), (aian(|)), tv - atan(^), TV + atan(^), 2tv — atan(^))) 

str(2) = (4, 16, (V32, V32, \/32, V32), (atan(|)), tt - atan(i), TV + atan(i), 2tv - atan(i))) 
str(3) = (4, 15, (V34, V34, \/34, \/34), (ata«(|)), tt - atan(f), 7r + atan(f), 2yr - atan($))) 
str(4) = (4, 12, (V40, ViO, VW, VW), (atan(|)), 7r - oian(4), 7r + atan(i), 2tv — atan{^))) 
str(5) = (4, 9, {VTS, vT8, yTS, vT8), (ntan(|)), TV — at<m(|), 7r + atan(|),2» — qtqn(j|))) 

2 6 40 str(l) = (4, 12, (V25, V25, V25, V25), (atan(i)), tv - atan('i), TV + atan(%), 2tv — otim(|))) 

str{2) = (4, 16, (-^32, s/32, v/32, n/32), (at hii(t)), 7t — atan(^j-J, 7r + atan( j), 27r — atan(^j-))) 
str(3) = (4, 15, (-/34, \/34, \/34, V^M), (aton(f)), tt - atan(^), TV + atore(f), 2yr - oton(|))) 
str(4) = (4, 20, (vTT, VTT, ViT, V2T), {atan(l)), n - atan(i), tv + atan(f), 2tt - atan(i))) 
str{5) = (4, 25, (\/50, \/50, \/50, \/50), (atonff )), tt - atan(^), tv + atan(i-), 2tv — <jtan(|))) 
str(6) — (4, 18, (V^, \f&5, V45), (atan(|)), 7r — atan(^), TV + atan(^), 2tv — atan(^))) 

3 9 40 str(l) = (4, 12, (\/25, v"25, \/25, V25), (atan(f)), tt - atan(4), tv + atan(4), 2tv - atan{4))) 

str{2) = (4, 16, (VS2, V32, V32, s/32), (ataraff )), 7r - atan(l), TV + atan(i), 2tv - atan(l))) 
str(3) = (4, 15, (v"34, V"34, VZi, %/Zi), (atara(I)), tt - atan($), tv + atan(|), 2n - atan(%))) 
str(4) = (4, 20, (v 7 !!, V^T, V^T), (atara(| )), it - atan(i), TV + atan(i), 2tv - otan(|))) 
str(5) = (4, 25, (V50, VS0, \/50, \/60), (atara(i)), tt - atim(4), TT + atan(4 ), 2ir — atan(|))) 
str{6) = (4, 12, (\/40, v'lo, v'lo, \/40), (atara(|)), it - oian(|), tt + ata»($), 2tv - atanA))) 
str{7) — (4, 18, (V45, V45, \/45, V45), (atan(4)), 7T — atan(^), TV + atan{|), 27r — atan(^))) 
str(8) = (4, 24, {V52, Vb2, s/52, Vb2), (atan(S)), it - otan(|), TV + atan(I), 2tt - atan(i))) 
3tr(9) = (4, 30, (n/61, n/61, n/61, n/61), (atan(^)), tt - atan(^), TV + atan(^) 2tv — atan(^))) 

4 20 100 str(l) = (4, 10, (V22.25, V22.25, V22.25, V22.25), (atan(^f)) , tv - aian(^i), tv + atan(^), 

27T — oton(^))) 

str(2) = (4, 8, (-v/20, ^20, ^20, V20), (at 

str(3) = (4, 15, (\/i54, \/34, \/34, \/34), (ata n(|)) , tt - otan(|), tt + aian(|), 2tv - otan(|))) 
str(4) = (4, 14, (V28.25, V28.25, V28.25, V28 25), (atan( g^)), 7r - oion( ), tt + atan( ^g), 
2tv — atan(^?- ))) 

str(5) = (4, 7^50, (V27.25, V27.25, V27.25, V27.25), (atan(ijS)), tv - atan(Jj£), tt + ata?i(Aj£), 
27T — atan(ij^))) 

str(6) = (4, 18, (V45, s/IS, VIS, \/4E), (atan(^)), tt - otan(|), TV + atan(|), 2tv - atan(|))) 
str(7) = (4, 12, (v^O, v^O, ^40, \/40), (atara(|)), tt - atan(S), tt + atan(§), 2tt - aian(|))) 
str(8) = (4, 18, (V45, V45, V45, v/45), (atan(|)), tt - atim(§ ), tt + atan(|), 2tt - oian(|))) 
str(9) = (4, 20, (v 7 !!, v 7 !!, Vil, \/4T), (atan(|)), tt - oian(|), tt + aian(|), 2tt - atan(i))) 
str(10) = (4, 5.25, (V14.50, V14.50, \/14.50, \/14.50), (atan( 53)), vr - aton(i^), tv + atan(^), 
2tt — atan(i^))) 

Str(U) = (4, 12, (v^5, V25, V25, V25), (otan( j)), 77 - atan(|), 7T + oton(|), 2ir - atan(^))) 
str(12) = (4, 6, (s/18.25, V18.25, n/18.25, n/18.25), (atan{^)),TV - atan(ii), 7r + atan(±£-), 
2tv — atan(ij^))) 

str(13) = (4, 15, (-v/34, ^34, ^34, \/34), (oten(|)), tt - atan(i), tv + atan(|), 2tt - atan(f ))) 
str(14) = (4, 20, (-v/IT, VE, v 7 ^!, \/4T), (atan(|)), tt - atan(|),ir + atan(|), 2w - atan(|))) 
str(15) = (4, 17.50, (V37.25, V37.25, VWM, V37.25), (atan(3jS)), tt - atan(^), 7r + atan(^), 
2tv - atani^f-))) 

str(16) = (4, 15, (V42.25, V42.25, V42. 25, V42.25), (atan{^-)), tv - atan(2^-),TV + atan(2Ji ), 
2tv — a(an(^))) 

str(17) = (4, 12, (ViO, v 7 ®, v 7 !?), V40), (atan(^)), TV - atan(^), TV + atan(^) , 2tv - atan(^))) 
str(lS) = (4, 20, (VJT, v 7 !!, VIT, V2T), (atan(|)), tt - atan(|), TT + atan(|), 2tt - atan(l))) 
str(19) = (4, 30, (n/61, \/61, V^T, v 7 !!), (atan(|)), tt - aton(|), 7T + atan{^) , 2tv - aton(|))) 
str(20) = (4, 9, (vT8, vT8, vT8, v 7 !^), (atan(j)), 7r - atan(^) , TV + aimi{|), 2tv - atan(^))) 



of the SA algorithm to terminate. When the problem size is increased to 20, the 
benefits of the SA algorithm begin to become apparent, as it out-performs the other 
two algorithms in terms of both solution quality and run time. In order to establish 
the significance of this, we now test all three methods on much larger instances. 



4.3. Large Rectangular Instances 

We designed instances with 40, 60, 80 and 100 rectangles. Space precludes a detailed 
description of these, but the full problem set is available from the corresponding 
author. As before, each method was run 50 times on each instance. Because of 
the computational cost incurred, we reduced CYCLE to 1000 for each algorithm. 
The results are depicted in Table El with an example solution for the 40 rectangle 
instance depicted in Figure d 

The SA method significantly out-performs the other two methods in terms of 
solution quality, but with an associated cost in terms of run time. However, as 
shown by the figures for standard deviation, SA offers a consistently high-quality 
solution method (at a price), whereas the other two algorithms offer solutions of 
more variable quality, but more quickly. 
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Table 4. Results for SA, CA-PSLS and PSO on four rectangular instances from the literature (CYCLE=3000) 



Instance 


Size 


Algorithm 


Tht>9+ 

' UCSL 


r 


(7 


t (s) 


1 


5 


SA 


12.776 


13.693 


0.62 


0.82 






CA-PSLS 


10.942 


11.704 


0.49 


4.31 






PSO 


11.046 


11.716 


0.49 


4.6 


59 


2 


6 


SA 


16.004 


17.377 


0.69 


1.66 






CA-PSLS 


14.686 


15.590 


0.67 


7.76 






PSO 


14.320 


15.349 


0.56 


8.28 


3 


9 


SA 


20.849 


22.328 


0.87 


6.84 






CA-PSLS 


18.157 


19.797 


1.03 


21.07 






PSO 


18.579 


19.205 


0.49 


22. 


84 


4 


20 


SA 


29.969 


31.680 


0.98 


92.01 






CA-PSLS 


27.927 


33.129 


5.48 


125.52 






PSO 


32.596 


34.426 


2.23 


138.00 


Table 5. Results for SA, CA-PSLS and PSO on large instances (CYCLE- 


= 1000) 




Instance 


Size 


Algorithm 


Hest 


r 


To 


I 00 


1 


40 


SA 


164.061 


174.586 


4.81 


263.24 






CA-PSLS 


179.508 


253.627 


60.42 


219.84 






PSO 


242.471 


276.939 


24.44 


197.03 


2 


60 


SA 


170.284 


187.312 


6.32 


905.75 






CA-PSLS 


184.984 


288.642 


124.26 


579.31 






PSO 


272.282 


317.739 


26.17 


451.72 


3 


80 


SA 


265.654 


281.087 


8. 


54 


2178.31 






CA-PSLS 


298.524 


544.421 


162.81 


1016.70 






PSO 


432.347 


490.862 


35.75 


813.59 


4 


100 


SA 


406.991 


423.087 


7. 


87 


4260.54 






CA-PSLS 


658.352 


880.537 


10* 


3.65 


1611.52 






PSO 


598.265 


688.785 


48.06 


1277.43 



5. Conclusions 

In this paper we describe a novel algorithm based on simulated annealing for the 
problem of packing weighted polygons inside a circular container. As well as being 
of significant theoretical interest, this problem has real significance in domains such 
as satellite design in the aerospace industry. Our algorithm consistently generates 
high-quality solutions that offer a significant improvement over those generated by 
other methods. However, this superiority comes with an associated computational 
overhead, so the choice of method should largely be driven by the anticipated appli- 
cation. Future work will involve improving the method's performance on problems 
containing nonconvex polygons, as well as its extension into three dimensions. 
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Figure 7. Best 40 rectangle solution generated by SA (r = 134.07) 
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